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TECHNICAL MEMORANDUM X- 64842 


SPECTRAL I RRADIANCE CURVE CALCULATIONS FOR 
ANY TYPE OF SOLAR ECLIPSE 

I. INTRODUCTION 


The spectral irradiance curve (SIC) gives the variation of the sun’s 
radiation of a given wavelength reaching a distant observer for various positions 
of the moon across the sun. It is simply the light curve for a given wavelength 
(A). The SIC is defined by the difference (l-o?) of the eclipse function (EF), 
cy , from unity. The EF gives the variation of the loss of solar radiation reach- 
ing the distant observer during the various phases of the eclipse. It is also 
referred to as the Russell Function in the theory of light curves for eclipsing 
binaries [l-3l and is defined in the List of Symbols. 

According to von Oppolzer's Canon of Eclipses 1 4 ] , an average of 
approximately 237 solar eclipses occur per century, of which 66 become total, 
77 become annular, and 10 are annular-and-total, and about 84 are partial at 
maximum. 

This report contains a simple procedure for calculating the eclipse 
function, and hence the SIC, for any type of solar eclipse and for radiation of 
any wavelength for which the solar limb-darkening obeys the limb-darkening 
cosine law: 


J=J (l-x + x cosy ) 
c 


( 1 ) 


where the symbols are defined in the List of Symbols. The case of occultation 
(partial/total) eclipses is presented in Section II and that of transit (partial 
annular) eclipses is discussed in Section III. 

The computation of the eclipse function or the SIC would be trivial but 
for the presence of the solar limb -darkening effect. The treatment here differs 
from that of Beard [5] which is based on the Method of Julius [6] and consists 
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essentially of dividing the sun into 12 somewhat arbitrary concentric zones and 
assigning' to each an average brightness value. The treatment in this report is 
based on the theory for the light curves of eclipsing binaries \ 1-3 1 , where it is 
assumed that the orbit is circular and that the stars are spherical and appear 
darkened at the limb according to the cosine law, equation (l). In the proce- 
dure here the binary system is replaced by the sun and the moon which are 
assumed to be spherical. Because the speed of the earth-moon system in its 
orbit around the sun is more than 30 times as great as the speed of the moon 
in its orbit around the earth, the relative orbit of the moon around the sun is 
always concave toward the sun and can be regarded as circular with respect 
to the sun throughout the entire interval of the eclipse. 

Also, the authors of this report consider only those solar radiations for 
which the limb-darkening law of equation ( l) holds for the sun. One can obtain 
the values of the limb-darkening coefficient (x) for these wavelengths from 
Reference 6. Then, the pertinent portions of Merrill's Tables [1-3] (herein- 
after called "Princeton 23") can easily be applied to the case of the solar 
eclipses. 

As an example of the computations involved, calculations have been 
made for the March 7, 1970, eclipse. One can use the procedure in a similar 
manner to compute the eclipse function and the SIC for any other occultation or 
transit eclipse. The March 7, 1970, eclipse is an example of an occultation 
eclipse which is total for a portion of the time, as seen from within a certain 
narrow track. 

Before proceeding to the description of the method, a brief explanation 
of the concepts of eclipsing binaries as applicable to the case of solar eclipses 
is given. 


A. Concepts and Relations Relevant to Solar Eclipses 

If 6 is the longitude in orbit from conjunction and i is the inclination 
of the orbit, then the apparent distance 6 between the centers of the two disks 
(of the sun and the moon) is given by (Fig. 1) 



cos 2 i + sin 2 i sin 2 0 


( 2 ) 
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Figure 1. Apparent geometry of eclipses. [The eclipse is an 
occulta tion if the eclipsed area ( shaded) is on the smaller 
disk. The eclipse is transit if the eclipsed area is 
on the larger disk. ] 

where f is simply a factor relating the distance units involved in the two sys- 
tems. The geometrical depth of the eclipse is then defined by 


p = (6 - r )/r , (3) 

g s 


where r and r are the respective angular radii of the larger and the 
g s 

smaller disks as seen by the observer. Depending on the distance of the moon 
from the observer (which varies about a mean distance, D ) , the moon’s 

disk may have its radius larger or smaller than that of the sun. When for any 
fixed observer the moon’s semidiameter (r ) is greater than that of the sun 

(r ) , the eclipse is herein called occupation; if and while p =£ -1, the eclipse 
sn 

is complete and is usually called total. When for any fixed observer the moon’s 

semidiameter (r ) is less than the sun’s (r ), the eclipse is herein called 
m sn 

transit; if and while p < _ 1 , the eclipse is complete and usually called annu- 
lar. When p > -1 , whether during occupation or transit, the eclipse is 
usually called partial. 



The definition of p applies equally, whether the larger or the smaller 
disk is in front: The eclipse will be absent, partial, or complete, according 
to whether p 5 l, 1 > p> -1, or p < -1 . Because the amount of light obscured 
during an annular phase (if there is one) is least at internal tangency (p - -l) 
and increases steadily with diminishing p to a maximum (loss) at concentric 
eclipse (p = -l/k) , if that occurs, it is more efficient to treat occultation 
eclipses as a whole and transit eclipses as a whole, by slightly different 
techniques. 

The other important quantity is 


k = r /r 
s g 


( 4 ) 


which is always less than unity unless r = r . From equations (3) and (4), 
one obtains [1-3] 


5 = r ( 1 + kp) 
g 


( 5 ) 


Only for a concentric eclipse (annular or total) , p c = -l/k (where k and p 

are dimensionless quantities), do whole values completely define the geometrical 
circumstances of a given phase. With the aid of geometrical parameters of 
the solar eclipse, the EF or the SIC is derived for the occultation eclipse in 
the following sections. 


B. The Eclipse Function and the Spectral I rradiance 
Curve Calculations 


The eclipse function, a , referred to as the Russell Function in the 
nomenclature of eclipsing binaries, is defined as the ratio of the light loss at 
any instant during the eclipse to the light loss at internal tangency. cv is a 
function of the limb-darkening coefficient, x, and the eclipse parameters k 
and p defined previously, i.e. , 
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a - a (x, k, p) 


( 6 ) 


One can invert relation (6) and obtain p as a function of x, k and ot , i.e,, 


P = p(x, k, a ) 


(7) 


during all partial and annular phases but not, of course, during actual total 
phases since ot is constant there. 

If a 0 is the value of a at mideclipse, then one can define the light 
loss at any time during the eclipses relative to that at mideclipse by 


n = ot /ot 0 


( 8 ) 


Then by definition. 


PO * P(X, k, Q! 0 ) 


(9) 


From equations (2) and (5) 


f(cos 2 i + sin 2 isin 2 0) 


r 2 ( 1 + kp) 2 
g 


( 10 ) 


For rr>idGclipse, 0 = 0 so that 

f 2 cos 2 i = r^ 2 (1 + kp 0 ) 2 

tr> 


( 11 ) 


Furthermore, if p = p(x, k, na 0 ) and 0(n) refers to the value of 6 when 
n 

the radiation loss from the sun’s disk is ot = na 0 , one can obtain the following 
important relation [ 1-3] : 
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( 12 ) 


sin 2 9 (n) 



(1 + kp n ) 2 - (1 + kp 0 ) 2 

= x(x, k, a 0 , n) 

(1 + kp x/ ) 2 - (1 + kp 0 ) 2 


or 


sin 2 9 (n) 


sin 2 9 (1/2) x(x, k, a 0 , n) . 


(13) 


These equations apply to all eclipses. It should be noted that n becomes unity 
at second contact (internal tangency) and that during totality the independent 
variable in x is p rather than n . 

The SIC, the light curve for a given X , is defined by the relation 
between 9 and (l - na; 0 ) for any eclipse. The SIC can be easily obtained by 
using Princeton 23 [1-3] along with equation ( 13 ) , as described in Sections II 
and III. 


II. THE OCCULTATION ECLIPSE 
A. The Total Eclipse 

For all occultati on-type eclipses, for all observers, k = r /r ^ 1. 

For an observer on the centerline of the path of totality, p 0 = - l/k ; for one on 
either edge of that path, p 0 = -1; for any observer elsewhere within that path, 

- l/k < p 0 < -1; and for one outside the path of totality but within the (very 
large) region where the eclipse is partial at maximum, -1 < p 0 < +1. For any 

observer within the path of totality a^ C = 1 and, in fact, o; OC = a° C = 1 through- 
out the time interval from second to third contact; for one outside that path but 

oc 

inside the large region mentioned, a 0 < 1 . The parameter p takes on the 
instantaneous values +1, -1, -1, +1, respectively, at first, second, third, and 
fourth contacts for an observer within the path of totality; second and third 
contacts are, of course, not observed outside that path, but the name "fourth 
contact” is by convention retained for the end of the phenomenon (Fig. l) . 
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Princeton 23 provides among other things x° C as a function of k 
and £V? C for a large number of selected values of x and n, and simple 

oc 

interpolating forms for finer screening in any one of the four variables k, or 0 , 
\ , n where warranted. 

An example of computation of 0(n) values is given below for the March 
7, 1970, centerline of totality for which p 0 = -l/k. From the Naval Observa- 
tory circular [7] , one obtains for an observer at the centerline of totality the 
value of Po = -1.040 , so that 


k° C = -l/p? C = 0.962 


(14) 


Then if the sun’s semidiameter 17] with respect to the earth's surface is 

r = 0. 004688 radians, that of the moon is 
sn 


r = r /k - 0. 0048773 radians , (15) 

m sn 


so that at first contact (p= +1, n=0.0) 

0( n = 0. 0) = r + r = 0. 00956451 radians . (16) 

' sn m 

From the x=0.6 section of Princeton 23 [ 1-3] , one can obtain 6 x° C (n) to 
use in equation (13) or, in a form more convenient for our present purpose, 
such as 


sin0(n) = sin 0(0.5) £ 6 x(n)J ^ 

to obtain 0{n) . The value of 0(0.5) is derived by substituting in equation 

(17) the value of 6 x° C (n= 0.0, k = 0.962, a 0 = 1, n= 0. 0) obtained from the 

next to the bottom row (x C ) on page 180 of Princeton 23, and that of 0(0. 0) 
from equation (16) : 



sinfl (0.5) = sin[0(O.O)]/[ O * 6 x(n)] ^ 


= sin (0. 0095641) /2. 30954 


=• 0.0041413 


( 18 ) 


Therefore, 0(0.5) = 0.0041413 , for n= 0.5 and (l -nor 0 ) = 0.5. Similarly, 
for n= 0. 1 and (l - no; 0 ) = 0.9): 


sin0(O.l) = sine (0.5) ^x(0.l)J ^ 


= 0.0041413 (1.83078) = 0.0075818 


0(0.1) = 0.0075818 


and so on. The results of the calculations for 9(n) and 0(l-na o ) are given 
in Table 1 and the SIC is plotted in Figure 2. The calculations have been made 
for Walterboro, S. C. (32 deg 54 min north, 80 deg 40 min west) , which is on 
the centerline of totality for the March 7, 1970, eclipse (see map in Reference 
7). The universal times [7] for the first and second contacts are 17 hr 4 min 
56 sec and 18 hr 24 min 67 sec, respectively. 

The time interval. At, elapsed after first contact can easily be obtained 
from 0's because of the uniform (circular) motion of the moon around the sun, 
over that interval. If the time interval between the first and second contacts 
is T 12 , and = r sn + r m [as in equation ( 17) ] , 



where 0 2 is the value at second contact. 
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TABLE 1. CALCULATIONS FOR 0(n) AND 0(l - na 0 ) [For Walterboro, 
South Carolina (32 deg 54 min north, 80 deg 40 min west) on the 
Centerline of Totality and x= 0.6 Corresponding to 
A = 550 nml 


Eclipse Phase 

n 

1-n 

e(n) 

(radians) 

At( n) 
( min) 

1st Contact 

0,0 

1.00 

0.009561 

0.0 


0.1 

0.90 

0.007 579 

16.780 


0.2 

0.80 

0. 006536 

25.617 


0.4 

0.60 

0.004883 

39.621 


0.5 

0.5 

0.004144 

45.880 


0.6 

0.4 

0.003427 

51.953 


0.8 

0.2 

0.001987 

64.149 


0.9 

0.1 

0.001210 

70.730 


0.95 

0.05 

0.0007794 

74.377 

2nd Contact 

1.00 

0.0 

0.0002125 

79.1784 


Note that, since 0 (n) is practically always so small during any solar 
eclipse (seen from a single point on the earth’s surface), ^(radians) is always 
equal to sin e(n) to the accuracy of the data. Its obvious implication would 

y 

be that a permanent set of tables of [ xl 2 for relevant values of x, k, n could 
very easily be constructed from the appropriate x's of Princeton 23. 

If H 0 (X, 0) represents the spectral irradiance ( watts/cm 2 -\) of the 
sun at instant t = 0 just prior to the first contact and H 0 (\, t) is that at time 
t after contact, the relative spectral irradiance is given by the ratio 
[ H 0 (A, t)/H 0 (A, 0) ] and is equal to (l - na 0 ). Hence the name SIC for the 
relation (l - nn 0 ) versus 9 or t. 
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RELATIVE SPECTRAL IRRADIANCE, 


1.0 


SPECTRAL IRRADIANCE CURVE 



Figure 2. Computed spectral irradiance curve for the limb-darkening 

coefficient, x= 0.6. 
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B. Occultation Eclipses Not Central at Maximum 


For an observer situated away from the centerline of totality, the 
relation p 0 = -l/k does not hold and one has to obtain the values of k and p 0 
either from Duncombe's circular [7l (if listed there directly) or by simple 
geometry. If the values are listed, then one proceeds as before; otherwise, 
one may proceed as follows. 

1 . To obtain k from geometry (March 7, 1970, eclipse) — For an 
observer at P (Fig. 3), let the moon be at a zenith distance, Z = 36 deg. The 
semidiameters of the moon and the sun given in the circular [7] are measured 
geocentrically and have the following values: R g = 16 min 6.8 sec, R m = 16 

min 32 sec, respectively. From the earth’s surface, the semidiameter of the 
moon is then given by the relation = R m D/ d where D is the instantaneous 

distance between the centers of the moon and the earth and is given by the 

relation D = D (0 /0), and d= D sin C/sin Z: 

m m 


r = R sinZ/sin(Z - M) 

(20) 

m m 


k = r /r = R sin(Z-M)/(R sin M) 
sms m 

(21) 


where 


sin M =• R„ sinZ/D 

ill 


( 22 ) 


For the March 7, 1970, eclipse, 0 = 60 min 39.4 sec and, therefore, D = 

359 464 km (224 665 mi) so that k = 0. 961. Thus, knowing the zenith distance, 
Z, of the moon at the time of the eclipse and knowing D, R g , R m and R E » 
one can determine k at any location. 

2. To obtain the value of Pq (March 7, 1970, eclipse) — In the eclips- 
ing binary notation, p 0 is the mideclipse value of p at the centerline. For 
observers on the northern and southern limits of the track of totality (Fig. 4), 
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c 


Figure 3. Geometry for deriving the instantaneous apparent 
diameter of the moon at a given place. 

p 0 = -1.000. The magnitude of eclipse values (m ) given by Duncombe [7] for 

J±i 

points outside the path of totality, i.e., for the cases of occultation eclipse 
partial at maximum, are related to the values of the eclipsing binary p 0 by 
the formula 


p 0 = (1-2 m £ ) 


(23) 
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Figure 4. Track of totality (or annularity) for the 
occultation (or transit) solar eclipse. 

If, however, the magnitude for a place is not given, one can obtain the 

value of p 0 by interpolating linearly with a distance from -1. 040 (at center- 
line) towards -1.000 (at the edge approximately 44 km away) for a place P 
inside the path, or extrapolating beyond -1.000 for a place not too far outside. 
The distances are measured perpendicular to the track. For places not too 
far away from the tracks, one can measure off perpendicular distances on the 
maps (in Duncombe's circular [7l) themselves. For example, consider the 
case of an observer at Tampa, Florida. Measuring the distances perpendicular 
to the track one obtains for Tampa: centerline (C.L. ) to the edge of the track, 
d = 23 mm, and centerline to Tampa, dp = 77 mm, so that for Tampa 


^ p ° 3 Tampa 




(lPol C.L. 


1 . 00 ) 


77 

= - 1.040 + — (0.040) = - 0.906 


(24) 


rs 



On the other hand, the magnitude of eclipse at Tampa, as given in the circular 

[7], is m = 0.953 so that 
E 

[pj = 1 - 2 (0.953) = - 0.906 , (25) 

F0 Tampa 


which is in agreement with the result in equation (24) . 

3 . To obtain the value of oJq needed — If P lies inside the path of 

totality, a 0 ° C = 1 whatever the particular values of x, k and p 0 . If P lies 
in the large outside region where the eclipse is partial at maximum, an 

oc 

assumption must be made as to the degree of darkening X ; then the needed 
is located in the appropriate a° C (k, p) table of Princeton 23. 

4. If the eclipse does not become total at any point of the earth's sur- 
face, then one resorts to the American Ephemeris and the Explanatory Supple- 
ment to the Ephemeris, or to local circumstances supplied directly on request 
by the Naval Observatory. 

5. The SIC — Given the angular radii of the sun and moon, the x, k, 

X a^ C and p° C , the determination of points (n, 0(n) or At), and therefore 
points on the SIC, proceeds as before. 

III. THE TRANSIT ECLIPSE 


In transit eclipses, the angular radius of the moon is less than that of 
the sun. Equations ( l) through ( 13) are still valid but one uses the tr tables 
of Princeton 23, instead of their occ counterparts, in the calculations, with 

= r /r ; while for occultation eclipses, = 1 at deepest phase for an 
m sn tr 

observer anywhere within the path of totality, in the transit case cnj >1 foi 

an observer within the path of annularity. 
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tx* fcl* 

At centrality of the annular case k = and p 0 = - !/ m E > where 

m is the magnitude of eclipse. In the large region of the earth for which the 

transit eclipse is partial at maximum, k is obtained by the geometry for the 

tl* 

phase (as for occultation eclipses) but Pq is obtained from either 

pj r = l _ 2 m /R or is read from a map as in the example of equation (25) . 

E 

All the procedures are essentially formally parallel throughout. Some 
relations appear different because 

tr _ light lost at a given phase of transit eclipse 


where r is the light that is or would be lost at internal tangency; r(x,k) is 
the tabulated in extemo in Princeton 23. 

During the annular phases of an eclipse observed from a point within 
the path of annularity, the amount of light lost at a given instant is strongly 
dependent on the degree of limb-darkening of the solar disk at the wavelength 
of observation. For this reason the small section of the light curve (or SIC) 
during annularity holds special possible interest. For more detailed study of 
this small section, Princeton 23 contains, in each x section a table headed 

X n 3nn . It is used by workers in the field of eclipsing binaries, as it would be 

fcl* 

here, in conjunction with the x tables for the same x > so it is placed with 
those tables. Its mode of use is obvious from its structure. 


IV. CONCLUSIONS 


A fairly direct method for obtaining the spectral irradiance curves for 
various assumed degrees of limb -darkening for any type of solar eclipse has 

1 / 

been indicated. As mentioned earlier, a permanent set of tables of I xl 2 
for relevant values of x, k, n could very easily be constructed from the 
appropriate x’s of Princeton 23, and the values could be stored in a computer 
data bank or on tapes for future use. 
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